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Abstract. The transverse momentum dependence of the anisotropic flow V2 for it, 
K, nucleon. A, S and O is studied for Au-|-Au coUisions at ^snn ~ 200 GeV within 
two independent string- hadron transport approaches (RQMD and UrQMD). Although 
both models reach only 60% of the absolute magnitude of the measured V2, they both 
predict the particle type dependence of V2, as observed by the RHIC experiments: 
V2 exhibits a hadron-mass hierarchy (HMH) in the low pT region and a number-of- 
constituent-quark (NCQ) dependence in the intermediate px region. The failure of the 
hadronic models to reproduce the absolute magnitude of the observed V2 indicates that 
transport calculations of heavy ion collisions at RHIC must incorporate interactions 
among quarks and gluons in the early, hot and dense phase. The presence of an NCQ 
scaling in the string-hadron model results suggests that the particle-type dependencies 
observed in heavy-ion collisions at intermediate pr might be related to the hadronic 
cross sections in vacuum rather than to the hadronization process itself. 
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1. Introduction 

The main goal of liigh energy nuclear collisions is the exploration of QCD-matter at 
high temperatures and densities. The search for a new form of matter consisting of 
deconfined quarks and gluons as constituents (usually called the Quark-Gluon-Plasma, 
QGP) has been a focal point of theoretical and experimental studies over the last decade 
pp. Measurements of the collective motion of hadrons produced in high-energy nuclear 
collisions have long been suggested as a valuable tool to gain information about the 
nature of the constituents and the equation of state of the system in the early stage of the 
reaction [3 El IH El El 13 El El lEl • Specifically, strange and multi-strange hadron elliptic 
and radial flow results seem to indicate that the observed collectivity originates from 
a partonic phasej JT]. Furthermore, elliptic and radial flow measurements for heavy- 
flavor hadrons like J / ^ and D mesons will test the hypothesis of early thermalization 
in these collisions P ^ n ^ iri]. 

At RHIC, measurements ^SlCSl of the elliptic flow V2 and the nuclear modification 
factor Raa,cp for identified particles have led to the conclusion that hadrons ought to be 
formed via the coalescence or recombination of massive quarks [T7 1 ITS l ITO l 120 1 l2T | l22 l l28 l 
EllESl- A cornerstone of this conclusion is the observed number-of-constituent-quark 
(NCQ-) scaling of the flow of baryons vs. mesons. Because this interpretation addresses 
key issues in high-energy nuclear collisions such as deconfinement and chiral symmetry 
restoration, it is of utmost importance to conduct a systematic study of other possible 
explanations for the observed particle type dependencies. 

For the present analysis we employ two independent hadron-string transport models 
RQMD(v2.4) and UrQMD(v2.2) [21 EH to study the effect of hadronic cross sections, 
kinematics etc. on the particle type dependence of V2. We present the V2 values of vr, 
K, p, A, S and Q for Au+Au colhsions at ^/s^^ = 200 GeV. We find that - although 
the hadronic models give only 60% of the magnitude of V2 - both models reproduce 
the gross features of the particle type dependence, including the mass ordering at 
Pt < I GeV/c and the number-of-constituent-quark dependence at > 1 GeV/c. 
Within both models, the NCQ dependence is related to the additive-quark-model for 
hadronic cross-sections. These findings imply that detailed comparisons between data 
and theory are necessary to disentangle hadronic and partonic effects at intermediate 
transverse momenta, 1.5 GeV/c < px < 5 GeV/c. For a transport model discussion 
of the quantitative accuracy and validity of the various methods to extract elliptic flow 
values from the data, the reader is referred to [H El Ej 

As one can observe from Fig. Qthe build-up of the elliptic flow in these models is 
tightly connected to the hadron-|-hadron scattering rate (top/bottom figures on the Ihs.). 
Initial valence (di-)quark scatterings before 1 fm/c do only yield a negligible amount of 
V2- The right panels in Fig. ^show the time evolution of the elliptic flow for mesons and 
baryons for low Pt (bottom) and high px (top). Already at early times, a separation 
between meson and baryon elliptic flow is visible, leading to the mass hierarchy at low 

I The term partonic is used to denote any kind of deconfined matter made of quarks and gluons. 
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Figure 1. UrQMD calculation for Au+Au interactions at -^snn — 200 GeV. Top 
left: Meson and baryon V2 versus time at midrapidity. Bottom left: Collision rates as a 
function of time for hadron+hadron (dashed line), hadron+quark/di-quark (full line) 
and quark/di-quark+quark/di-quark (dotted line). Bottom right: Meson and baryon 
V2 versus time at midrapidity and pt < 1 GeV/c. Top right: Meson and baryon V2 
versus time at midrapidity and pT > 2 GeV/c. 



transverse momenta and to the stronger elliptic flow of baryons (compared to mesons) 
at high Pt- 

2. Model Results 

Within the framework of the hadronic transport approach, a typical heavy ion collision 
proceeds schematically in three stages, i.e. the pre-hadronic (strings and constituent 
(di-) quarks) stage, the hadronic pre-equilibrium stage, the evolution towards hadronic 
kinetic equilibrium and freeze-out. The pre-hadronic stage involves the initial excitation 
and fragmentation of color strings and ropes. At the highest RHIC energy, this stage 
lasts for about 0.5-1.5 fm/c and the effective transverse pressure/EOS is rather soft. 
During the late hadronic stage, the hadronic system approaches local kinetic equilibrium 
followed by an approach to free-streaming, where the system escapes equilibrium due to 
dilution of the hadronic gas: the mean free path of the hadrons exceeds the finite size of 
the system [221 HZl I2H] , the free streaming hadrons decay and feed down to the lightest 
species. 



Anisotropic flow at RHIC 



4 



0.08 - 



0.06 



0.04 - 



0.02 









Ai 1 


+ Au ai 






* STAR 


:, 


j 


★ 


★ 


★ 






RQMD(v2.4) 










★ 




• default 
O no res. 






A 


A 


A 




★ 


UrQMD(v2.2) 
A default 


"a 


A 


• 


• 


• 


A 
• 


A 
• 


★ 

A 
• 

A 




--■0-- 




--0 


0" 


.Q 


— a 





100 200 300 400 

Number of participants Np^^j 



Figure 2. Charged hadron V2 versus the number of participating nucleons in Au+Au 
interactions at y^snn = 200 GeV. Experimental data OD] from the 4-particle 
cumulant method is shown as stars. RQMD results are depicted by filled circles (full 
calculation) and open circles (without rescattering). The UiQMD calculations are 
shown as open triangles. 



Figure |21 presents the model results on the centrality dependence of the charged 
hadron f 2- values along with measurements from the STAR collaboration IHO] . Both 
hadronic transport models (UrQMD v2.2 and RQMD v2.4) reach about 60% of the 
measured V2 values only, although the centrality dependencies are very similar to the 
data. There is a small variance between the two models, which we consider as an 
estimate of the systematic errors in such model calculations. Although the V2 values 
from the hadronic transport model also depend on the formation-time of hadrons from 
strings [7j, the failure of both hadronic transport models to describe quantitatively 
the magnitude of V2 is a strong indication that there are interactions amongst pre- 
hadronic constituents (partons) present in nature (but not in the hadron-string models 
discussed here), which are responsible for the large V2 values observed in the experiments 
inn IS21 inni- when rescattering between the hadrons is turned off (full circles), V2 
vanishes completely, because repulsive vector interactions are not included into the 
present simulations [3 El IZ| • 

Let us now analyze the temporal structure of the elliptic flow's development. 
Figure El shows the f 2 values as a function of the freeze-out time for pions in minimum 
biased Au+Au interactions at ^snn = 200 GeV. The different symbols denote different 
transverse momentum intervals decreasing from top to bottom. One clearly observes 
a strong correlation between freeze-out time and elliptic flow: particles that decouple 
earlier have a larger V2 value than those that freeze-out later. The lower px particles 
tend to freeze-out later and their V2 continues to evolve late in the evolution (a similar 
correlation has been observed earlier in low energy collisions |3^) The decrease in the 
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Figure 3. Elliptic flow {v2) of pions from RQMD versus their freeze-out time for 
several pt windows for minimum bias Au+Au collisions at y^snn = 200 GeV. The 
dashed-line represents the pr integrated distribution. 



V2 towards later times is related to the reduction of the coordinate-space anisotropy 
with time - i.e. the system has become more spherical than it was at earlier times. The 
f 2 values of higher p^-pions saturate sooner and tend to reflect the earlier stage of the 
collision more strongly. Thus, within the model dynamics, the final V2 is mostly driven 
by the early stages of the reaction; the V2 values at high px are closer to the initial/early 
V2 values than the V2 values at lower px- 

3. Particle Type Dependence 

How much of the observed NCQ-scaling features can be reproduced by the hadronic 
models? In both dynamical approaches, finite (vacuum) cross sections are used to model 
the strong interactions in the hadron-string cascade. Unlike the simplistic Cooper-Frye 
freeze-out treatment in most hydrodynamic calculations, the transition from strongly 
interacting matter to free-streaming is determined here by the interplay of the local 
particle density and the energy dependent cross section of the individual hadrons. It is 
well known that a proper treatment of the gradual freeze-out is crucial for the finally 
observed hadron distributions. It was pointed out that the hydrodynamical results on 
flow depend strongly on the proper kinetic treatment of the freeze-out process and can 
not be approximated by isotherms [23 lEHl IHZj ■ 

However, the major shortcoming of the present hadron-string approach is the 
lack of the early partonic interactions which are important for the early dynamics in 
ultra-relativistic heavy ion collisions El El 13 IHH]- In order to take care of both 
partonic and hadronic interactions in high-energy nuclear collisions, a combination 
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Figure 4. RQMD results of tt, fC, p, and A V2 from Au+Au collisions at ^snn = 
200 GeV. (a) Minimum bias collisions: At about pr ^ 1-2 GeV/c, baryon and meson 
V2 are crossing each other; (b) 40-80%; (c) 20-40%; (d) 0-20%. 



of the hydrodynamic model for the early stage dynamics (the "perfect" fluid stage) 
with a hadronic transport model for the later stage plus freeze-out has been proposed 
jnni Sni HH in ESI EH- Fig. m shows the colHsion centrality dependence of the pr- 
dependent V2 values for vr, p, and A. Both, the hydrodynamic behavior (in the low 
Pt region) and a hadron-type dependence (in the intermediate pt region) are clearly 
predicted in all centrality bins. This "crossing and subsequent splitting" between meson- 
and baryon elliptic flow as well as the breakdown of the hydrodynamical mass scaling 
at high transverse momenta was first predicted within the UrQMD model [7] and has 
later been observed in the experimental data. It is important to note that the more 
recent explanations of this effect (the suggested "number-of-constituent-quark" scahng) 
is not a unique feature of the "quark recombination/coalescence" assumption: hadronic 
interactions alone have quantitatively (at the correct pr-values) predicted this hadron 
type dependence. 

Let us explore the pr-dependence of the event anisotropy parameters in detail. 
Figure El shows the calculated unsealed and scaled f 2 values of various hadrons versus 
the unsealed and scaled transverse momenta, pt, of the various hadrons. On the left 
hand side (FigureEfa)) one can see that at lower transverse momenta, px < 1-5 GeV/c, 
the heavier hadrons exhibit smaller V2 values than the lighter hadrons: Hadron transport 
theory predicts mass ordering. The S and fl V2 values from the UrQMD calculations 
are also included. They are the lowest of all V2 values for pt < 2 GeV/c. Such mass 
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Figure 5. (a) Hadron V2 from minimum bias Au+Au collisions at ^snn = 200 GeV; 
(b) Scaled hadron V2 results are shown. The Uq refers to the number of constituent 
quarks. Symbols represent results from the UrQMD (v2.2) model for various hadron 
species. At low px/Nq < 0.5 GeV/c, tt does not follow the scaling perhaps caused by 
the resonance decay 1X4}. In higher pr region, K meson seems to fall off the scaling 
curve due to the comparatively small hadronic cross sections in the model. 



ordering is exactly what is observed in the experimental data ^3] and is in accord with 
hydrodynamic calculations [45j. Hence, hadronic interactions, which do take place at 
later stages of the collisions, also do contribute to the observed collective motion. 

At higher px values, this mass dependence gives way to the f2(pT)-dependence 
on the hadron type (i.e. meson or baryon). Here, it is interesting to note that the 
r2-baryons seemingly acquire a significant amount of V2 in the model calculations. In 
addition, there is also clear, but small difference for kaons and pions in V2 values at 
Pt ^ 1-5 GeV/c. This particle type dependence, rather than the otherwise dominating 
particle mass dependence, is also observed in the data (TH] . It is important to note that 
the meson has a mass that is very close to the mass of the baryons p and A, and, 
indeed, recent experimental results on the 0's V2 values are similar to other mesons |46j . 
However, in the hadronic transport model, about 2/3 of the 0-mesons are formed via 
ii'--S'-coalescence, which is not necessarily the dominant process in heavy ion collisions 
|47j . Therefore, the V2 values of meson are not shown in Figure El It should also be 
noted that in the high px region, px > 2.5 GeV/c, all V2 values start to decrease. This 
indicates that the system is deviating from an ideal hydrodynamic behavior. This trend 
is best seen in the right, "scaled", plot in Fig. Such a drop has been observed in the 
data. 

The test of the NCQ scaling hypothesis is shown in Figure Efb), which depicts the 
scaled hadron values, V2/nq. The scaling factor is the number of constituent quarks 
(NCQ) in accord with the coalescence approach |2SII31I2S]- For mesons and baryons. 
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riq = 2 and ^^=3, respectively. The NCQ-scaling is clearly observed in both RQMD (not 
shown here) and UrQMD model calculations except for the pions. This surprising result 
and its implications for the frequently invoked recombination/ coalescence hypothesis will 
be discussed in the last section. 

However, one should note that there are experimentally distinguishable differences 
(when differential Vt elliptic flow data with good statistics becomes available) between 
"real" NCQ scaling from coalescence and the approximate NCQ scaling due to the cross 
section hierarchy discussed here. In the coalescence approach one expects identical 
elliptic flow of all baryons irrespective of their type. In contrast, the present AQM 
scaling picture, leads to an observable ordering of the elliptic flow at fixed pt according 
to the strangeness content. This can indeed be observed in Fig. |31 where an ordering 
of the elliptic flow values can be observed: f2(N) > f2(A) > f2(S) > V2{^)- 

4. Discussion and Conclusions 

The particle and energy density is highest at the center of the created fireball in 
relativistic nuclear collisions - initially, there is an angular dependent matter density 
gradient. The repulsive interactions among the constituents will therefore push matter 
to move outwards. In this way, the collective flow develops in nucleus-nucleus collisions 
[H 121 ini HI m E] • We would like to stress that flow means matter- and energy flow. It 
is independent of the type of particles, either partons or hadrons, or different kinds of 
hadrons. Hence, by studying the collective motion of the produced hadrons one can, in 
principle, extract the information of early collision dynamics jU El • In general, one 
expects that the final elliptic flow. 



where S denotes the hyper-surface where hadrons are emitted, will depend on a, i.e. 
the interaction cross section, which, in principle, depends on the particle type, cm angle 
and relative momenta. The specific particle density depends on the collision time t, 
location, and momentum. For short mean free paths, the transverse flow is intimately 
related to the pressure, which in turn depends on the density and temperature of the 
matter under study [H El El IH El El IZ| • Indeed, the frequent rescatterings among the 
hadrons can lead to hydrodynamic-like mass ordering in the low p^ region. 

At the higher transverse momenta, pt > 1.5 GeV/c, the particles escape quickly 
from the system to low density, in effect leading to early freeze-out and lack of 
development of the hydrodynamics and the details of the interaction cross-sections are 
most important. As the cross sections depend on the particle type, for mesons or baryons 
to first approximation given by the constituent quark model we do expect roughly 
a 2:3 scaling of the meson-to-baryon elliptic flow from transport calculations. 

The hadronic models underpredicted the strength of V2 at RHIC, because early 
partonic interactions (except from quark coalescence during the string break-up) are 
not included in the model. The early stage with highest density and smallest mean free 




(1) 
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paths is "missing". This shortcoming of the hadronic models clearly demonstrates the 
need for the early, dense partonic interactions in heavy ion collisions at RHIC. 

5. Summary 

In summary, the hadronic transport models UrQMD and RQMD have been employed to 
study the elliptic flow of hadrons in Au+Au collisions at the highest RHIC energy. We 
have analyzed the V2 values as a function of collision centrality, transverse momenta and 
collision time for various meson- and baryon species, including multi-strange baryons. 
Both hadron-string transport models fail by 40% to exhaust the absolute value of 
V2i probably due to the lack of (partonic) interactions in the initial, hot and dense 
stage. However, because of the hadronic re-interactions, the hadron mass hierarchy is 
qualitatively well reproduced in the low pt region. Rescattering is the key that leads 
to the quasi- ideal hydrodynamic appearance in V2{pt)- Also in the intermediate pt 
region, the hadron type dependence (number-of-constituent-quark scaling) is predicted 
by both hadronic transport models. Here, this dependence is due to the hadronic cross 
sections which do roughly scale with the number of constituent quarks, in accord with 
the additive quark model. This finding challenges the recent interpretation of NCQ 
scaling as a unique deconfinement signature. Thus, further tests of the deconfinement- 
plus-recombination hypothesis are necessary with high precision V2 measurements of 
resonance hadrons like K*, p, A, A*, and S*. 
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